Classical dimers on the triangular lattice 
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We study the classical hard-core dimer model on the triangular lattice. Following Kasteleyn's 
fundamental theorem on planar graphs, this problem is soluble using Pfafhans. This model is par- 
ticularly interesting for, unlike the dimer problems on the bipartite square and hexagonal lattices, its 
correlations are short ranged with a correlation length of less than one lattice constant. We compute 
the dimer-dimer and monomer-monomer correlators, and find that the model is deconfining: the 
monomer-monomer correlator falls off exponentially to a constant value sin(7r/12)/\/3 w .1494 . . . , 
only slightly below the nearest-neighbor value of 1/6. We also consider the anisotropic triangular 
lattice model in which the square lattice is perturbed by diagonal bonds of one orientation and small 
fugacity. We show that the model becomes non-critical immediately and that this perturbation is 
equivalent to adding a mass term to each of two Majorana fermions that are present in the long 
wavelength limit of the square lattice problem. 



I. INTRODUCTION 

The study of classical dimer models - the statistical 
mechanics of hardcore dimer coverings of graphs - has 
a venerable history. These models have been of inter- 
est as direct representations of the physics, e.g. diatomic 
molecules on a lattice but even more because of their 
equivalence to various other statistical mechanical prob- 
lems; for example, the two-dimensional Ising model can 
be reformulated as a dimer model on a special latticei^ 
The reduction to dimer form is advantageous in that a 
wide class of dimer models, those on planar graphs with 
independent fugacities on bonds, are soluble by Pfaffians 
following a theorem of Kasteleyni 3 . This technique has 
been applied to compute a variety of correlation func- 
tions, in the dimer model on the square lattice^ and in 
the Ising model. 5 

The interest in dimer models has received new and 
very different impetus following the discovery of high- 
temperature superconductivity. Following Anderson's 
proposal 6 that the superconducting state evolves out of 
a liquid of singlet (valence) bonds, Rokhsar and Kivel- 
son proposed a quantum dimer modeli to describe this 
so-called short-range resonating valence bond (RVB) 
physics. The Hilbert space of their model consists of 
all pairings of the spins on the square lattice into singlet 
bonds; these can simply be labeled by hardcore dimer 
coverings. The quantum dynamics provides off-diagonal 
matrix elements between these coverings. This idea can 
be generalized to other lattices. 

The properties of the classical dimer model enter the 
solution of such quantum dimer models in two ways. 
First, trivially but usefully, the infinite temperature stat- 
ics of the quantum problem are given precisely by the 



classical problem with equal fugacities. Second, non- 
trivially and even more usefully, quantum dimer mod- 
els generically exhibit an "RK" point - which general- 
izes Rokhsar and Kivelson's construction on the square 
lattice - where the quantum wavefunction is an equal 
amplitude superposition of dimer coverings which is the 
paradigmatic RVB form. Static, dimer diagonal, ground 
state correlations at the RK point are then again given 
by the correlations of the classical dimer model£& For 
the square lattice, the known results on the classical 
problem 4 showed that the quantum dimer model was 
critical at the RK point, which thus turned out to be an 
isolated critical point between two solid phases, rather 
than a representative of an RVB phase. Similar behav- 
ior occurs on the honeycomb lattice, where the classical 
model is equivalent to the five- vertex model on the square 
latticed 

In recent work, two of the present authors have shown 
that the triangular lattice leads to a a different outcome 9 
with an RVB phase, including the RK point, character- 
ized by liquid correlations. It was also noted that the 
model exhibits "topological order" in the sense of Wen 11 
(ground state degeneracies on closed surfaces in the ab- 
sence of symmetry breaking) characteristic of an Ising 
gauge theoryAiS As part of this work it was necessary 
to solve the isotropic classical dimer model on the trian- 
gular lattice and show that the dimer correlations were 
short-ranged. As this does not appear to have been done 
previously, with the exception of results for the thermo- 
dynamic limit entropy^&iiiiSii&ii we have carried out a 
fuller analysis of the problem in this paper. 

Two aspects of this expanded analysis are notewor- 
thy. First, we show that monomers are deconfined on 
the triangular lattice. A monomer is a site which does 



2 



not have a dimer touching it, and the classical monomer- 
monomer correlator is the ratio of the number of config- 
urations with the two monomers to the number without 
them. Monomers are deconfincd if this ratio approaches 
a constant, nonzero value as the distance between the two 
increases. In the quantum dimer model, two monomers 
have the interpretation of two "test spinons" obtained 
by breaking a valence bond apart into its constituent 
spins. Hence our result proves that spinons are decon- 
fined at high temperatures in the triangular lattice quan- 
tum dimer model^AS Second, we study the interpolation 
between the critical square lattice and the non-critical 
triangular lattice by tuning the fugacity for one species 
of bond. The diverging correlation length near the for- 
mer allows a continuum limit to be taken and yields a 
theory of two Majorana fermions with the monomer op- 
erator being identified with a linear combination of the 
spin operators in the two sectors. 

We begin by computing the partition function for gen- 
eral fugacities and show that the model undergoes phase 
transitions only when it reduces to the square lattice 
model. Correspondingly, it exhibits long-ranged dimer 
correlations in the square lattice limit. Next we derive 
asymptotic forms for the Green function and dimer cor- 
relations which show that the correlation length is less 
than a lattice constant at the isotropic point. We ex- 
tend the analysis to obtain the correlation length every- 
where along the interpolation to the square lattice. In 
Section IV we turn to computing the free energy of two 
monomers in a background of dimers. We show that it 
falls off exponentially to a constant, proving that spinons 
are deconfined in the quantum dimer model at infinite 
temperature. In Section V we show that the continuum 
formulation of the critical square lattice problem is a the- 
ory of two Majorana fermions. The addition of the re- 
maining bonds of the triangular lattice is equivalent to 
adding a mass term for each fermion and hence a tran- 
sition of the Ising universality class. We end with some 
concluding remarks and two technical appendices. 



II. MODEL AND PARTITION FUNCTION 

A dimer is a bond connecting two nearest neighbors 
on a lattice. We study the close-packed model with hard 
cores, where an allowed dimer configuration has the prop- 
erty that each site of the lattice is paired with exactly 
one of its nearest neighbors, such a pair being denoted 
by a dimer placed on the link between the two sites. In 
the simplest form of the model, each dimer has the same 
fugacity; as the number of dimers is the same in all con- 
figurations, the correlations of the dimer model are thus 
given as the equally-weighted average over all possible 
dimer configurations. In the following, we will include 
unequal fugacities, so that the average to be taken then 
includes non-trivial weighting factors. 

The close-packed hard-core dimer model can solved on 
any planar lattice by using Pfaffian techniques. (A pla- 



nar graph contains no overlapping links.) These tech- 
niques were introduced in the early 1960s by Kasteleyn, 1 
and Fisher and coworkers A22i2i The result is simple to 
describe: on a planar graph one places arrows on the 
links so that each plaquette is "clockwise odd" , that is 
to say that the product of the orientations of the arrows 
around any even-length elementary plaquette traversed 
clockwise is odd. The antisymmetric matrix My is then 
defined with My = 1 if an arrow points from point i to 
point j, Mij = — 1 if the arrow points from j to i, and 
My = if i and j are not nearest neighbors. Kaste- 
leyn's theorem states that for any planar graph M can 
be found, and that the number of dimer coverings (the 
partition function) is given by the Pfaffian of the matrix 
M: 

Z = ±Pf[M]. 

The ± sign is chosen to make Z positive; henceforth we 
will omit this sign. This result is exceptionally useful 
because the Pfaffian is the square root of the determinant: 

Pf[M] = (det[M]) 1/2 . 

For regular lattices, this determinant can easily be com- 
puted by using Fourier transformation. This result 
has been reformulated in the more modern language of 
fermionic path integrals in Ref. ^[ One places a Grass- 
mann variable "0i on each site i of the lattice, and defines 
the action S = J2i<j ^iji^i^j ■ Then a basic result of 
Grassmannian integrals gives 

Z = J [Vtp]exp(S). 

This form allows one easily to define different but equiv- 
alent M by rescaling the fermions. 

For the triangular lattice, an appropriate choice of ar- 
rows is displayed in Fig. ^ To make contact with the 
work on the square lattice, we however choose a different 
convention. First, we deform the triangular lattice into a 
topologically (but not symmetry) equivalent square lat- 
tice. Second, we multiply the Grassmann variables on 
every other row by alternating factors of ±i so that the 
values of My are given as in Fig. ^ where the weight is 
to be understood to be for a convention of arrows point- 
ing to the right/upwards. Finally, in order to be able to 
interpolate between square and triangular lattice, we al- 
low the fugacity of the diagonal bonds to be a variable t 
instead of 1. For completeness, we give fugacities v and u 
to horizontal and vertical bonds, respectively. Note that 
we must double the unit cell to contain two sites; in our 
convention, it is doubled in the vertical direction. Each 
site i is labeled by the location of its unit cell, x, and its 
location in the basis, a, while the unit vectors between 
cells are x and y. The two fermions in the unit cell lo- 
cated at x are thus denoted as ij) a ,x with a = 1,2. In 
this form, the action is 

S =jEE M x%*y (1) 
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FIG. 1: Choice of M for the triangular lattice, (a) A 'clock- 
wise odd' sign convention, (b) By multiplying the Grassmann 
variables on the sites by the phases indicated (italics) , one ob- 
tains the new values for the Mij on the bonds. Bonds are in 
addition multiplied by a fugacity. Note that in both cases, 
the unit cell is doubled. 



where M™^ = —M@£. With periodic boundary condi- 
tions (or far away from the boundaries) , we have Mg£ = 
M Q/3 (x - y), where 

Af 12 (0) = M 21 (y) = iu, 
M n (x) = M 22 (x) = v, 
M 21 (x + y) = -M 12 (x)=it 

Since the action is quadratic in terms of the fermions, 
the model can be solved by Fourier transformation. Our 
convention for Fourier transforms is 



The action is then 



S 



o X] M k' 3 V'k,QV'-k,/3 



(2) 



where the two-by-two matrix is 

M = (2ivsmk x g(k) 
k 9* s ) 2ivsink x 
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TABLE I: Number of dimer coverings, Z, of triangular lattices 
of with L x x 2L y sites, with periodic and open boundary 
conditions. For periodic boundary conditions, only those sizes 
(L x > 3 and L y > 2) are given in which any pair of sites is 
linked by at most one bond. 



site, S, of the dimer coverings on a TV-site lattice is then 



S=llnPf[M] = i]Tln|A(k)| 



(4) 



with 2 



7(k) = i[u — te 



ik. r 



ue 



Finding the determinant and hence the Pfaffian is now 
simple, because in k space the action is expressed in 
terms of the four-by-four blocks 



(3) 



on the diagonal. Note that this matrix is antisymmetric 
as it must be, because <?*(k) = — g{— k). The entropy per 




where 

A(k) = det[M k ] = det[M_ k ] 

= -4v 2 sin 2 k x - 4u 2 sin 2 (fey/2) 
-At 2 cos 2 (fez + ky/2) 

For a system with periodic boundary conditions, one 
must sum over four different sectors, according to the 
winding of the dimer configuration^ In fermionic lan- 
guage, this corresponds to evaluating four Pfaffians, with 
symmetric and antisymmetric boundary conditions, and 
combining them as Z = (— Pf pp + Pf ap + Pf pa +Pf aa )/2, 
where the subscripts denote the (anti)periodic boundary 
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conditions for the fermions in the two directions. In a 
Fourier representation, this corresponds to the usual dif- 
ferent choices of allowed wavevectors. We have listed the 
number of dimcr coverings thus obtained for the isotropic 
case in Table U along with the results for a system with 
open boundary conditions. 

In the thermodynamic limit, the sum over k in Q 
turns into an integral. For the isotropic triangular lattice 
(t = it = v = 1), doing the integral numerically yields 

5 = 0.4286.... 

We note S has been obtained previously by several 
authorsiiiiiSiiSiii in different contexts, including the 
kagome Heisenberg magneU& and, implicitly, the fully 
frustrated Ising model on the hexagonal lattice^ where 
the ground state entropy per spin Sffhim = S/2 
is related to the dimcr model entropy via a duality 
mapping^ 

This model simplifies in several limiting cases. If two 
fugacities vanish, the model reduces to decoupled chains, 
with order along the chains and disorder relative to one 
another. More interesting is the case where one of the 
fugacities vanishes, the square-lattice dimer model. This 
model is critical, with algebraic decay of correlators^^ 
It is straightforward to show that the entropy S(t,v,u) 
is non-analytic in t at t = (and likewise for u and v). 
We will discuss this in more detail in section V. 



III. THE GREEN FUNCTION 

We now turn our attention to the correlations of the 
triangular dimer model. We utilize the techniques de- 
veloped for the square lattice in Ref. El and expressed 
in terms of Grassmann variables in Ref. [lfj Correlation 
functions can be straightforwardly expressed in terms of 
the Grassmann variables (see e.g. Ref. I24)) . Since the 
action is quadratic, using Wick's theorem expresses all 
correlators in terms of the Green function 



'Vip] ipiipj exp(S). 



For example, the probability P^ d \j) of finding a dimer 
on a bond in the (j) direction at site is given by 



P {d) 0) = \(M ] ) 



(5) 



In the original matrix language, the Green function ele- 
ments make up the inverse of the matrix M, namely 

= (M- 1 ) ji = -(M- 1 ) ij . 

In this section we explicitly determine the Green func- 
tion and its asymptotic behavior. In the next we use 
these results to determine the dimer-dimer and monomer- 
monomer correlators. 

As shown in the last section, in Fourier space the ac- 
tion is written in terms of the four-by-four blocks J2J). 



These blocks can be easily inverted to give the two-point 
functions in Fourier space. Expressed in terms of the 
functions <?(k) and A(k), we have 



(V'l.kV'l.-k) = (V'2,kV'2,-k) 



2i sin(fe x ) 
A(k) 



(V^k^.-k) = 
(V>2,kV>l,-k) = 



g*(k) 
A(k) 

g(k) 
A(k) 



To determine the Green functions in real space, we need 
to invert the Fourier transform. This yields 

(^,x#,y) = y"dk e - 4k -( X -^(^ k ^_ k ) (6) 

where we define 



dk 



Air 2 



2tt 



dk.. 



2tt 



dk*/, 



For a finite number of sites, this integral is replaced with 
the usual sum. 

To simplify the resulting expressions, we specialize 
slightly in the subsequent analysis: 

1. We work in the thermodynamic limit so we have 
integrals over k. 

2. We set u = v = 1. Thus the fugacity t gives a 
way of interpolating between the square and the 
isotropic triangular lattices, with t = giving the 
former and t = 1 the latter. 

3. We study correlators where the fermions are in the 
same or adjacent rows. 

It is also convenient to absorb some factors of ±i by defin- 
ing the Green functions Q s and R s as 

(V>l,xV>l,x+ S x) = (V>2,xV>2,x+ S x) = (-1) (S+1)/2 Q S 
(V>l,x^a,x+«*) = -(V^xVl.x-sx) =1(-1) 1{S+1)/2] R S 

where [a] is the greatest integer less than a. Relatively 
simple expressions for Q s and R s are obtained by shifting 
k x ~ > k x + 7r/2 and k y — > 2k y + tt. For odd s 



<ik w(k x , k y ) cos(k x ) cos(sk x ) (7) 
R s — t J dk w(k x , k y ) cos(k x + k y ) cos(sk x + k y ) (8) 

while for even s 

Qs = (9) 
R s = J dk w(k x , k v ) cos(ky) cos(sk x + k v ) (10) 

where 



w(k x ,k y ) 



1 1 

2 cos 2 (k x ) + cos 2 (k y ) + t 2 cos 2 (k x + k y ) 
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These integrands are all invariant under the interchange 
of k x and k y , under k — ► — k, and under shifts of ir. 

One of the two integrals in each of the Q s and R s can 
be done immediately by residue. We let a = k x + k y and 
b = k x — k y so that u> _1 = 2(l + cos(a) cos(6) + i 2 cos 2 (a)). 
The integral over b can be deformed so that the contour 
runs around the pole at 



cos(fe) 



1 + t 2 cos 2 (a) 
cos(a) 



This is done by changing the integration contour running 
from — 7r to 7r along the real axis to one which first runs 
from — 7r to — 7r + zoo, then to +7r + ioo, and finally back 
down to 7r. The two contours running parallel to the 
imaginary axis cancel as the function is periodic under 
b — ► b + 27r, and the contribution from the contour run- 
ning parallel to the real axis vanishes when the integrand 
vanishes exponentially as Im(fo) — > oo. We thus pick up 
only the contribution of the pole of the integrand. For 
example, we have 



db 



1 



1 



2vr 1 + rcos(6) +t 2 r 2 ^1 + (2t 2 - l)r 2 + i 4 r 4 ' 



The resulting expressions for Q s and R s are simpler 
on the isotropic triangular lattice t = 1 because of the 
additional symmetry. In particular, here w(k x , k y ) is in- 
variant under the transformation k x — > ~{k x + k y ); this 
ensures that the Green functions are invariant under 60- 
degree rotations of the lattice. For example, this trans- 
formation allows us to change the numerator of Q s from 
cos(sk x ) cos(k x ) to cos(sa) cos(a). This gives us 



Qs(t = l) 



and for odd s 



1 



da- 



cos(sa) cos(a) 



2tt J y/l + cos 2 (a) + cos 4 (a) 



R s (t 



while for even s 



R s (t = l) = 



1 

4-7T 



cos(sa)e i6 ° + cos((s - l)a) 
yjl + cos 2 (a) + cos 4 (a) 

cos((s — l)a)e lb ° + cos(sa) 
o \J\ + cos 2 (a) + cos 4 (a) 



da 



da 



(11) 



(12) 



(13) 



where 



ibo 



1 + cos 2 (a) — \Jl + cos 2 (a) + cos 4 (a) 
cos(a) 



A. Asymptotic long-distance behavior 

In this subsection, we study the long-distance behavior 
of the Green functions Q s and R s . We derive an expres- 
sion for the correlation length of the fermions at any t, 
and find that it indeed vanishes only when t = 0. 
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FIG. 2: The integration contours (dashed) used for the 
asymptotic evaluation of the Green function when t > 1/2. 
The second is equivalent to the first if the integrand is multi- 
plied by 2. For t < 1/2, the branch cuts have Re(a) = 0. 



Unlike the similar case of the square lattice, the inte- 
grals for the Green function cannot be evaluated asymp- 
totically by partial integration. This method only works 
for algebraic correlations, whereas in our case, the cor- 
relations decay exponentially. One can nonetheless ob- 
tain the asymptotics by deforming the integration con- 
tour. Instead of there being poles in the complex a-plane, 
there are now square-root branch cuts. In the region 
|Re(a)| < 7r/2, there is one branch cut in the lower-half- 
plane and one in the upper. For t < 1/2, the branch cut 
has Re(a) = 0, while for t > 1/2, it is parallel to the 
Re(a) axis. The general formula for the location of the 
branch points a± is 



cos(a±) 



1 ± Vl - 4t 2 
2t 2 



(14) 



For the same reason we could deform the 6-integral 
around the pole, we can now deform the a-integral around 
the branch cut, as depicted in first part of Fig. [5] To 
evaluate the integral by saddle point, we want a contour 
where the e zsa term in the integrand is decreasing ex- 
ponentially. It is thus easiest to use the fact that with 
a square-root branch cut, the integral around one side 
of the branch is equal to the integral around the other. 
Thus the contour can be deformed yet again, so that it 
goes from the two branch points to ioo, as long as we 
multiply the result by 2. This contour, denoted by C, is 
depicted in the second part of Fig. |21 

With this new contour, asymptotic evaluation of the 
integral is easy, because the integrands are sharply 
peaked at the branch points. For example, we have 



,(*=!) 



1 



da- 



O i(s+l)o _i_ „i(s-l)a 



2tt Jc y/l + cos 2 (a) +cos 4 (a) 



The exponentials rapidly decrease as Im(a) is increased 
along the contour, and moreover, the rest of the integrand 
has a square- root divergence at the branch points. It is 
thus an asymptotically exact approximation to substitute 



a/1 + cos 2 (a) + cos 4 (a) w C^/a — a± into the integrand. 
This gives 



Q.(t=l)«/(|* + l|) + /(|a-l|) (15) 



where 



/(*>) = Re 



Ctt 



da- 



11 + 



V° = a + 



We have used the fact that the integral from a_ to infin- 
ity is the complex conjugate of the integral from a+ to 
infinity. After shifting a — ► a + a+, the integral is easily 
done, yielding 



f(s) = Re 



1 



VCs 



Finally, we break a + into its real and imaginary parts: 
a + = a + i/3. For t = 1, cos(a + ) = exp(2z7r/3), so that 
C = 3 3/4 2 exp(-nr/4), a = 1.1960 ... and 



/3=iln 



y/2 + 3 1 / 4 
V2-3V4 



= .83144 . 



Plugging this in gives 
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(16) 



(17) 



The asymptotic expression for R s can now easily be found 
by noticing that at the branch point a = a+, e = — 1. 
This means that 



R s (t = i)n(-iy[f(\s\)-f(\ s -i\)] 



(18) 



We have therefore shown that on the isotropic trian- 
gular lattice, (3 is the inverse correlation length along 
the rows. It is only about one lattice spacing, so the 
exponential decay of the Green functions is quite rapid. 
Comparing the asymptotic form to the exact values de- 
rived in the next subsection (see Fig. [3]for a plot of Q s ), 
one finds that the agreement is excellent — within a few 
percent — even for small values of s. At occasional val- 
ues of s (e.g. s = 15), however, the agreement is only 
within a factor of 2 or so. The reason is that the oscil- 
lating factor in l|17l) occasionally becomes very close to 
zero near s an odd integer, so that the terms we have 
neglected above can become larger than the one we kept. 
This however, does not change our results for the corre- 
lation length, because the neglected terms have the same 
exponential dependence on s (but a power law exponent 
different from -1/2 in Ea. lT7|) . 

For arbitrary t, the computations are similar but the 
equations look fairly gruesome. Finding asymptotic ex- 
pressions for Q s and R s is more complicated because 
the transformation k x — ► — (k x + k y ) no longer leaves 
w{k Xl ky) invariant. However, we can easily extract the 
correlation length for Green functions in the x + y direc- 
tion (the direction along the links with ^-dependent fu- 
gacity) . The reason is that the expressions for the Green 
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FIG. 3: The curve is the asymptotic expression for \Q S \ when 
t = 1 l|15[l . The points are the exact values at odd s; Q s 
vanishes for even s. 



functions in this direction end up very similar to H12I13JI . 
The key fact is that there are no terms in the integrand 
involving e isb ° , only e ib ° . Like before, the integrand is 
peaked around the branch points given in (jl4J) . Thus the 
fermion correlation length £ x +y is given by 



Im 



arccos 



1-y/l- 4t 2 



(19) 



We plot this correlation length in Fig. Note that the 
kink at t = 1/2 in Fig. 0] is evidence for a level crossing. 
As t — > 0, the correlation length diverges as 1/t: 

1 



t -> o 



while £ = 1//3 = 1.2027 ... at t = 1. In this expression, 
we have taken the length of a diagonal bond as unit dis- 
tance. Translating this back into square-lattice language 
(with the length unit given by a lattice constant) removes 
the factor of l/y/2. The diverging correlation length in 
the square lattice limit will be used to take a continuum 
limit in Section V below. 



B. Short-distance behavior on the isotropic lattice 

In this subsection we show how to derive exact expres- 
sions for the Green functions on the isotropic triangular 
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FIG. 4: The correlation length as a function of the dimer 
fugacity: t = is the square lattice, t = 1 the isotropic trian- 
gular lattice. 



lattice, t = 1. We start by noting that a few of these 
Green functions can be evaluated immediately. In the 
isotropic case, the probability of a dimer being on a given 
link is 1/6. This is of course related to the Green function 
for neighboring fermions by (J5J so 

Qi = Ri = Ro = 7- 
6 

Ironically, this fact is not easy to extract from the explicit 
integrals (|1 111 2ll3fl . It does follow from noting that the 
transformation k x — — (k j: + k y ) in (|7I8I1UI) means that 
Qi = R x = R , and also that Q 1 + R 1 + R = 1/2. 

By generalizations of this argument, it is possible to 
derive recursion relations relating Q s and R s . For the 
square lattice they are discussed in Ref. |23- For the 
isotropic triangular lattice, the relations relating all the 
R s to the Q s can be derived using simple trigonometric 
identities. One finds here that 



R 



2j+2 



R 



2j+l 



R 



•2.1 



R-2j-i + Q 



2j+l 



and 



2R 



2j+l 



-2R 



2j 



Q2j+i — Qij-i + Sjo- 



The Kronecker delta in the latter identity arises from 
the integral J dk cos(jx) = Sjq. 



The recursion relations 



Qs 



1/6 

-1/2 + 1/6^ 



Uqtt 



8 U + 5/6 - 5/6 
-56 U + f + 7/3^| 
216 % -f + 5/3^| 

1288 tj , 1397 _ 385 V3 
5 U ° "r 6 6 U tt 



00 



29232 U, 



o 



_ 3443 
6 

2885 _ 
2 



2509 



6 C/ ir 
30970 

Jh2L 
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numerical value 



0. 166666666666666666667 
-0.024550581458226763328 
0.0022132632430383043346 
0.0000693679214938995513 
-0.000076203880846465563 
0.0000170664326460443169 
-0.000002068500172729645 
0.0000000086255295387778 



TABLE II: The Green function Q s explicitly 



involving only Q s are trickier to obtain; we discuss them 
in the Appendix. 

The recursion relations mean that Q s and R s for any 
s can be expressed in terms of Qi, R2 and Q3. These in 
turn can be evaluated in terms of elementary functions; 
in particular the combination 



U 



£(7/6) 
2r(2/3)V^ 



(20) 



As detailed in Appendix A, along with Q± = 1/6 we have 



R2 
Qs 



-2U 
1 1 Vs 



6 U n 



Using these facts, we can evaluate the Green functions 
to arbitrarily good accuracy numerically by using e.g. 
Maple. Although it is easy to iterate the recursion rela- 
tions on the computer, we have not succeeded in finding 
a closed- form expression for arbitrary Q s and R s . We 
have collected some of them in Table [n] Note that even 
though the coefficients are increasing exponentially with 
s, the three terms in each cancel almost perfectly to give 
the exponential falloff in s. 



IV. CORRELATORS 

A. The dimer-dimer correlation function 

The dimer-dimer correlation function is easily ex- 
pressed in terms of the Green functions. The operator 
Vh.xV^.x creates a dimer at site x pointing in the y direc- 
tion. The probability that there are two parallel dimers 
in the same row a distance s apart is therefore 

P {dd) {s) = (^l.xV^x^l.x+sxV^x+sx) 

By Wick's theorem, this decomposes into 

P( dd \s) = (Ro) 2 - {Qs) 2 - RsR-s 
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This therefore decays exponentially to its asymptotic 
value with half the correlation length of the Green func- 
tion. For the isotropic lattice, one can easily plug in the 
asymptotic forms (|15I18|I derived in the last section, and 
find the correlation length 1/(2/3) « 0.6014, less than one 
lattice spacing. 



B. The monomer-monomer correlation function 

As noted in the introduction, the asymptotic behavior 
of the monomer-monomer correlation function is quite 
interesting physically. A zero value means the model is 
confining, while a nonzero value means that the model is 
in a deconfincd phase. We will show in this subsection 
that on the isotropic triangular lattice, the latter is true. 

A monomer placed on a lattice site forbids a dimer 
from being placed on any of the links connected to the 
site. The monomer-monomer correlator P (mm) (<7,r) is 
the ratio of the number of configurations with monomers 
at sites q and r to the number of configurations without 
the monomers. Thus computing a monomer-monomer 
correlator follows from the partition function of the lat- 
tice with the two sites (and all the links connected to 
these sites) deleted. Since such a lattice is planar, Kaste- 
leyn's construction is still applicable. The one compli- 
cation is that we must ensure that on the lattice with 
deleted sites, the number of arrows is still clockwise odd. 
With the assignment in Figure 1, one sees immediately 
that the number of arrows around a deleted site is even. 
Thus this assignment must be modified by reversing one 
of the arrows in the plaquette around the deleted site. 
Reversing the arrow on this link then ruins the clockwise- 
odd assignment around the other plaquette this link bor- 
ders. Thus we must reverse one of the other arrows on 
this other plaquette. This in turn ruins another assign- 
ment, and so on. We thus must build up a string of 
reversed arrows, which must stretch from one monomer 
to the other. As long as the arrows are chosen to make all 
plaquettes clockwise odd, the monomer correlator is inde- 
pendent of the choice of path of the string. For monomers 
in adjacent rows, this construction is illustrated in Fig. 
[SJ links with the thick lines are part of the string. 

This is very much like constructing the spin-spin cor- 
relation function in the Ising model in terms of fcrmionic 
variables. In fact, on the square lattice, the monomer- 
monomer correlator was showni*2£ to have the same long- 
distance behavior (falloff as the square root of the dis- 
tance) as the spin-spin correlation function in two de- 
coupled Ising models. The precise relation between the 
two correlators for the square lattice was given in Ref. 
EjH In section V, we will extend this relation away from 
the square lattice. 

We define as the monomer-monomer correlation func- 
tion p( mm ) as the ratio of the number of configurations 
with the monomers to the number without. This can be 
written in terms of the fermions. Deleting a site merely 
corresponds to inserting a fermion at that site. Changing 




3' 4' 



FIG. 5: The string in the monomer-monomer correlator with 
p = 4. 

the sign of the arrow corresponds to changing the sign of 
the term of the fermion in the action. In our theory with 
action Ylj<k Mjkipjipk, for a monomers at sites q and r, 

P( m "%r) = \ (^ q Y[(l-2M ab M b )A)\ (21) 

where the product is over all the links connecting sites q 
and r which have reversed arrows. We specialize to the 
case of one monomer being at the origin and the other 
in an adjacent row p lattice spacings apart, as illustrated 
in Fig. [3] Labeling the fermions in the lower row 1 . . .p 
and those in the upper row as p + 1 . . . 2p gives 

P (mm) (p) = |(Vi(l-2^ p+1 ^ 2 )(l + 2^2V P +2) 

...(1-2^-1^)^)1 (22) 

Since the theory is free in terms of the fermions, one 
can use Wick's theorem to express this correlator as a 
product of Green functions. Evaluating it by brute force, 
however, is too difficult for all but low values of p, because 
the number of contractions grows exponentially. For a 
two-site separation, 

P (mm) (2) = R 2 (l - Ro) - 2R\ - 2Ql 

Using the results of the previous section for the Green 
functions, one has for t = 1 

P (mm) (2) = 4f/ /3 - 1/9 = .14657672599 . . . 

To proceed further, one must in revert to the Pfaffian 
methods developed for the square latticed Our calcula- 
tion generalizes these methods to the triangular lattice. 
The monomer-monomer correlator (|22H can be expressed 
in Pfaffian language as 

p(mm) = Pf(M (mm) ) 

Pf(M) 

where M^ mm ^ is the Kasteleyn matrix with the monomer 
sites removed, and the signs reversed on the links of the 
string connecting the monomers. The advantage of this 
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3 (mm) 



0. 166666666666666666666667 
0. 146576725991984081282220 
0.150263558036976118604604 
0. 149354528957809826733084 
0.149441157286260650652237 
0. 149430105031348455969637 
0.149429091105386654359851 
0. 149429312376742966925414 
0. 149429243123620452470446 
0. 149429245483558265741555 
0. 149429245691495770851350 
0. 149429245319212097747918 
0.149429245371125328892038 
0. 149429245361581959940523 
0. 149429245361256645973930 



TABLE III: The monomer-monomer correlation function for 
the isotropic triangular lattice 



formulation is that one can manipulate a matrix without 
changing the determinant. For example, if one adds any 
row of a matrix to any other row, the determinant is left 
unchanged. Thus the strategy is to find explicitly the 
matrix M _1 M( mm ) describing the monomer correlator, 
and manipulate it to make it more tractable^ On the 
square lattice, it is simple to show that 

det(Af- 1 M( mm )) = det(T) 2 , 

where T is a pxp matrix. After a variety of manipulations 
described in appendix B, one finds the same behavior on 
the triangular lattice. 

Deferring the details to appendix B, our result is that 
the monomer-monomer correlator is 



p(«w0(p) = I det(K + Q). 



(23) 



where the entries of the pxp matrices Q and 1Z are 

Ha = (-i^-^Rj-i+i+f-^eii-j) (24) 
Qa = i(-i) lu+t)/2] Q P+ i^-j (25) 

where 8(x) — 1 for x > and for x < 0. Unlike 
the case of the square lattice, the matrix 1Z + Q is not 
Toeplitz (the entries in a Toeplitz matrix depend only 
on i—j). For Toeplitz matrices, Szego's theorem provides 
a simple way of finding the asymptotic behavior of the 
determinant as the matrices get large. The theorem is 
not applicable here, because while IZij depends only on 
i—j, Qij depends on i + j. 

To evaluate this correlator for t = 1, we can plug the 
Green functions derived in the previous section into the 
expression for p( mm ) (|23|l . Expressing them in terms of 



Uq as before, the final form of the determinant simpli- 
fies substantially (although it is still fairly horrible). For 
example, for p = 3 one finds 



p( mm )(3) = 



V3 + (2tt - 3V3)C/ - 144tt[/ 3 + 864irU$ 



27ttU, 



o 



We have collected some numerical values in Table IIIII 
The correlator falls off exponentially to the value 

p(™™)(oo) = 0.14942924536134225401731517482693. . . 

To all 30+ digits, this number is given by 26 



p(™«)(oo) 



sin(7r/12) 
y/3 



Since this number is fairly simple (and involves sin(7r/12), 
which has already appeared in our analysis), we expect 
that this expression is exact and there exists a general- 
ization of Szego's theorem applicable to the determinant 

(J23J. 

The non- vanishing of the monomer one-point function 
(= \/P( m ™)(oo)) means that the monomers are decon- 
fined in this phase. Numerically, one can check that the 
correlation length here is the same as for the dimer-dimer 
correlation, namely 1/(2/3), or about .6 of a lattice spac- 
ing. 



C. The vison-vison correlation function 

In the quantum dimcr model, a 7r-flux vortex or 
"vison"^Zi 2 £ involves a semi-infinite string on the dual 
lattice. A particular dimer configuration c acquires a 
factor of (— l) String ( c ) : where JV s t r i n g(c) is the number of 
dimers intersecting that string. For two visons separated 
by x the string can be taken to run from one vison to the 
other and hence no flux is visible at infinity. At the RK 
point, the calculation of the static two vison correlator 
reduces to the evaluation of ~^2 C {— l) A '=t«ng<o,=)( c ) which is 
now a purely classical problem. We will henceforth refer 
to this sum and its generalizations to unequal fugacities 
as the two- vison correlator even though the latter do not 
correspond to any known quantum problem. We should 
also note that the choice of string, which needs to be 
made for each separation, lead to an ambiguity of a sign 
at each separation. 

The two-vison correlator was studied recently in 
Ref. |2^ for a Kagome spin problem equivalent to a gen- 
eralized dimer model on the triangular lattice requiring 
three dimers per site. There are some special features of 
this problem, arising from the "particle-hole" symmetry 
of occupying three dimers out of a possible six at each 
site and among them is the presence of two species of 
visons (for details the reader should consult the original 
work). In our problem, that of one dimer per site on 
the triangular lattice, life is simpler. Here the two-vison 
correlator can be seen, most simply by duality arguments 
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described in Appendix A of Ref.0, to be equal in magni- 
tude (we have already commented on the sign ambiguity) 
to the spin-spin correlation function averaged over the 
ground state manifold of the fully frustrated Ising model 
on the honeycomb lattice. The latter is known, and de- 
cays exponentially with a correlation length computed in 
Ref Ji. A direct computation of the vison correlator is in 
progress 313 and should obviate the need for the detour we 
have taken in this section. 



V. PERTURBING ABOUT THE SQUARE 
LATTICE: FIELD THEORY 

As we noted earlier, the general triangular lattice 
dimer model is critical in the square lattice limit. This 
suggests that we can gain some insight into its proper- 
ties by considering a continuum limit when the diago- 
nal fugacity t is small. We should immediately enter one 
caveat: our exact computation shows that there is a level 
crossing in the excited state spectrum of the transfer ma- 
trix (the kink in Fig. en route from the square lattice 
to the isotropic triangular lattice so there are details of 
the latter that the continuum theory will not reproduce 
very well. 

Returning to the fermion action and setting t = 0, it 
turns out that in the absence of diagonal bonds we no 
longer need to double the unit cell as was necessary for 
in Eq. [3J the fermion operators ipi used in this section 
thus refer to the sites x of the square lattice. 

In momentum space, the resulting action is given by 



So = 7T sink x — 2 sin fcj,)^^- 



(26) 



where we have now reverted to the standard Euclidean 
field theory convention of including an explicit minus sign 
with the action in writing down the Grassmann integral. 
We see that the fermion dispersion has nodes at k = 
(0,0), (0,7r), (tt, 0) and (it, ir). Linearizing near these 
points and defining 



xi 

~i 
Xk 

xl 

~2 

Ak 



il'-k 

^•^(x+yj+k 



(27) 



The leading asymptotic behavior of the lattice Green 
function can be recovered by inverting Eq. 1271 

V>(x) = X \x) + i-l) v x\x) 

+ i (-ir+^+ i x 2 (x^(-ir+ i x 2 (x) 

and using the continuum correlators 

1 1 



<X a (x) X a (0)) = 
<X a (x) X a (0)) = 



4-7T x + iy 
1 1 



4ir x — iy 
for a = 1,2. This procedure yields 

47r L_ x + iy x — iy J 

which exhibits the two sublattice structure that charac- 
terizes the square lattice problem. From this form we can 
recover the asymptotics of the dimer correlations estab- 
lished in Ref. 0- For example, the connected correlation 
function for two horizontal dimers displaced by x is given 
by the Grassmann expression, 

c xx (k) = (^(x)^(x + x)^(0)V(x)) - 

MxMx + *)W(0M&)> . 

By Wick's theorem, c xx = -G 2 (x) + G(x + x)G(x - x) 
and using our asymptotic form for G(x) we find 



Cxx (x) 



(— i) x +y 

4tt 2 



x + iy 



x — iy 



, (30) 



where we have ignored constants in the denominators 
that generate subleading corrections in the generic case. 
For the special case of two dimers in the same column, 
separated by an even distance, this expression vanishes. 
However, even in that case one can recover the correct 
subleading form by reinstating the constants, as can be 
readily verified. 

Next, we consider adding in the diagonal bonds with 
fugacity t and find that near the nodes of Sq it adds 



S p = it\ (-l)^(x)^(x + x + y) 



it ^2 exp(-ifc x ) cxp(-i(k y - 7r))V'k'0-k+7ry 



we obtain the long wavelength form (d = ±d x - id y , Expanding this near the nodes of the dispersion yields 
d = \d x + idy) 



So 



d 2 x (xW + X X dt + X 2 d X 2 + X 2 d X 2 ) (28) 



S p = 2it I d 2 x(-xY +X 2 X 2 ) 



which is a theory of two Majorana fermions, both holo- 
morphic (x^X 2 ) an d anti-holomorphic (x^X 2 ) fields. 
The labels 1 and 2 here are not the same as those used to 
label the two fermions in the unit cell; now the variables 
x and y run over all sites. 



(31) 



i.e. two mass terms of opposite signs. This means that the 
model is invariant under the standard Kramers- Wannicr 
duality of the Ising model, which sends t — > — t. Here 
this just interchanges the two fermions. It is easy to 
check that our Pfaffian analysis gives results independent 
of the sign of t. Changing t — > — t leaves invariant the 
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Green functions Q s , and R s for even s, while it flips the 
sign of R s for odd s. This leaves all monomer and dimcr 
correlators invariant. 

It follows that perturbing away from the square lat- 
tice results in a non-analyticity in the thermodynamics, 
with behavior in the d = 2 Ising universality class (with 
a doubling of the degrees of freedom). Specifically, On- 
sager's results for the Ising model give the dimer entropy 
density to be S(t) — 5(0) oc — t 2 \ogt 2 , so that the model 
is indeed non-analytic at t = 0. For t near zero, the 
correlation functions are now characterized by decay on 
the scale of £ ~ 1/t, in agreement with our earlier, exact 
lattice computation (JT^J. 

It is interesting to ask where the order/disorder fields 
(a J n) of the Ising operator algebra appear in the lat- 
tice problem. The two known candidates here are the 
monomer and vison fields, both of which exhibit an in- 
verse square root decay on the square lattic o 4 i 31 . In 
the Ising model at the critical point, order and disorder 
two-point functions decay as |x| 1 ^ 4 , so the monomer 
and vison fields must be bilinear in the order and dis- 
order fields, as proven in Ref. To narrow down 
the correspondence, we invoke four further constraints: 
(a) at t = 0, the monomer-monomer correlator is non- 
zero only when the monomers are on different sublat- 
tices (b) when t ^ 0, monomer correlators decay to a 
constant (c) the vison-vison correlation function does 
not exhibit any sublattice structure and decays exponen- 
tially to zero when t (d) under duality t — > —t, 
a and /i change places in each of Majorana sectors. Be- 
cause of the opposite signs in (|31fl . one of the sectors 
is in the high-temperature phase, while the other is in 
a low-temperature phase. In an Ising low-temperature 
phase (ct(x)ct(O)) goes to a constant at large |x|, while in 
a high-temperature phase (/x(x)/x(0)) goes to a constant. 
Thus for t ^ only one of the four bilinears decays to 
a constant, either a 1 /! 2 or a 2 /! 1 , where the superscripts 
refer to the two Majorana sectors. With these we can 
identify the monomer operators on the two sublattices 
with a 1 ^ 2 ± ia 2 /! 1 ; the relative ±i ensures that the cor- 
relator vanishes at t = when the monomers are on the 
same sublattice. In the absence of further constraints, 
the vison operator can be identified with either of the 
combinations ± [i 2 - A fourth dimension 1/4 op- 
erator still needs to be identified to complete this set of 
identifications. 

Finally, we should note that in the standard Ising 
model on the square lattice, equivalent to a dimer model 
on Fisher's lattice, things work somewhat differently, as 
we shall discuss elsewhere^ In that problem the two 
Ising phases correspond to confincd/dcconfined phases 
of the dimer problem. 



VI. CONCLUDING REMARKS 



fined monomers. This result contrasts with the criti- 
cal correlations found on the bipartite square and hon- 
eycomb lattices. The more complicated, non-bipartite, 
Fisher lattice also exhibits a deconfined phase along with 
a confincmcnt-dcconfinemcnt transition. It would clearly 
be interesting to see if other, as yet unstudied, examples 
follow this classification. In this context it is worth noting 
that the two critical dimer models admit height represen- 
tations on account of their bipartite character — a feature 
often associated with criticality in two dimensions. It is 
also worth remarking that this correlation between criti- 
cality and bipartedness connects, via the quantum dimcr 
model, to the Fradkin-Shenkcr theorem in lattice gauge 
theory; 33 Finally, for reasons noted in the introduction, 
any further work along these lines would also be immedi- 
ately useful in advancing our understanding of quantum 
magnetism. 



Note added 

After we submitted this paper, a preprint by 
A.Ioselevich, D.A.Ivanov, M.V.Feigelman appeared 
(cond-mat/0206451), who studied the triangular dimer 
correlations in the absence of monomers. Beyond the 
results reported here, they point out that the dimer cor- 
relation length depends on direction. 
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APPENDIX A: RECURSION RELATIONS 

To obtain the recursion relations for the Green func- 
tions efficiently, it is helpful to do changes several of 
variables. First of all, we define z = cos(2a). Since 
s is always odd, one can rewrite 2 cos(so) cos(a) = 
cos((s+l)a)+cos((s — l)a)) in terms of powers of cos(2a). 
Namely, we use Chebyshev polynomials Tj (y) , which are 
defined by the relation Tj(co$(a)) = cos(ja). For exam- 
pie, Tx{y) = y, T 2 (y) = 2y 2 - 1, T 3 (y) = 4y 3 - 3y, and 
so on; a closed-form expression can be found in 3 *. Using 
this gives 



We have shown that the classical model on a triangu- 
lar lattice is generically in a liquid phase with decon- 
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This means that any Q s can be expressed as a sum of 
the functions 

U k = ^- [ dl 



27T J_ X ^/(l-Z 2 )(7 + 4z + Z 2 ) 



APPENDIX B: DERIVATION OF THE MATRIX 
FOR THE MONOMER CORRELATOR 

The monomer- monomer correlator ()22l) is written in 
terms of matrices as 



.(*+ 


l)/2. 


For example, 






p( mm )(p) = Pf {M~ l M {mm) ). 


'l = 


Ux + 






We define 




'3 = 


2U 2 








E = M {mm) - M 


'3 - 


AU 3 


- 3U 2 + 2Ui - 






'4 = 


8L/ 4 - 


f 4U a - 8U 2 - 


aui + u 


Since by 


definition the Green's function 



One also has 



R 2 =2U X . 



The reason we write things in terms of the Ui is that 
it is easy to derive recursion relations for them. This is 
useful because there are simple explicit expressions for 
the first few U^. Denoting the denominator of (|A1|) as 
Y, the recursion relations follow from the fact thali^ 



1 

Y 



mz m - 1 Y' 1 



1 m d(Y 2 ) 



2" dz 

is a total derivative. For our case, this yields the identity 

{k - l)U k + 2(2k - 3)C/ fe _i + 6(fc - 2)C/ fe _2 

-2(2* - 5)C/ fc _3 - 7(k - 3)(7 fe _4 = 0. 

Even though the expressions for X s and Uk look un- 
wieldy, they can be made much nicer by a simple change 
of variable. Namely, defining z = (£ — 4)/(£ + 2) yields 



1 f° /f-4 



-in 



1 



1 



(A2) 



These turn out to be related to the complete elliptic in- 
tegral K(k) with k = sin(7r/12), an integral whose prop- 
erties were originally investigated by Legendre^ We can 
evaluate Uq in terms of the j3 function as 



Uo = 



£(7/6) 

2r(2/3)V^ 
.19326587782732139. . 



(A3) 



Since Q\ = 1/6, this yields immediately that 

£/i = 1/6 - E7b. (A4) 

This then gives us our first non-trivial Green function 
explicitly, because R 2 = 2U X = -.05319842232.... To 
utilize the recursion relations given above, we also need 
U 2 . Using various properties of the complete elliptic 
integrals^ we have shown that remarkably enough, U 2 
is simply related to the reciprocal of Uo , namely 



Uo = Un 



1 



4ttV3J7o 



(A5) 



(M 1 ) J - fc = -(ipjip k ), we have 

p(mm) = (det(/ + G£;)) 1/2 

where I is the identity. As seen from 12211 , there are three 
types of entries in E: 

1. We need to remove the links connected to the 
monomer sites 1 or 2p. Thus Ejk = —Mjk if ei- 
ther j or k is either 1 or 2p. 

2. We need to change the sign on the links on the 
string. Thus Ejk — —2Mjk if the link (jk) is part 
of the string. 

3. We need to account for the ipiip 2p in l|22l) . This is 
done by E 2p ,\ = ~Ei :2p = i. 

We now specialize to the case of monomers in adjacent 
rows a distance p apart. Here the matrix E has non- 
zero entries in 2(p + 4) columns and rows. In E a b, the 
indices r and s run over 1 . . . p for the sites in the lower 
row, and p + 1 . . . 2p for the upper row. The remaining 
8 columns and rows correspond to the sites surrounding 
the two monomers which are not part of the string; we 
denote those adjacent to the left monomer as 1', 2', 3', 4', 
and those adjacent to the right monomer as 5', 6', 7', 8'. 

Since we evaluated all the Green functions in section 
III for t = 1, at this point we could just plug in all the 
numbers, multiply the matrices, evaluate the determi- 
nant, and take the square root to get the monomer cor- 
relators. However, the determinant can be simplified a 
great deal, and is quite elegant. In fact, we can rewrite 
the correlator as a determinant of a matrix half the size. 

Let us first consider the first column of I + GE. 
The non-zero entries of Ebi are when b = 2p or b = 
2,p+ 1, 1', 2', 3', 4' (the latter sites surrounding the site 
1). Note that ^2 c G ac M c i = S a \ when c is summed over 
the sites around 1. Then 



GabEbl — iG a , 2 p — G ac Mcl — G a ^ 2p 



Sal 



Thus 



(7 + GE) al = iG a , 2p . 



Similarly, (I + GE) a , 2p = —iG a i- We now consider the 
columns coming from the sites adjacent to site 1 but not 
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part of the string. For these sites d = 1', 2', 3', 4', we have 
E a b $0,1- Thus in the A:th column of I + GE we have 
GabEbc' — G a \E\ c i. However, we can essentially remove 
this column without affecting the determinant. If we add 
the 2pth column times iE\ c i to the c'th column, the c'th 
column of (/ + GE) is just 5 ac i. Thus the only non-zero 
entry in the c'th column is on the diagonal, and is 1 . The 
determinant with these rows and columns removed is thus 
identical to those with these columns present. The rows 
and columns 5', 6', 7', 8' are removed in a similar manner. 

We denote the modified 2p x 2p matrix (with the same 
Pfaffian) as M. We can continue with such column ma- 
nipulations to simplify the matrix further. Consider the 
pth column. We have 

(I + GE) ap = S ap + 2iG a ,2 P -i — itG a £p- 

The last piece can be removed without changing the de- 
terminant by adding t times the first column of M to 
this column. This gives 

Map = Sap + 2iG a .2p-\- 

Now consider the p — 1th column: 

2iG a 2 P -2 — 2itG a ,2p-i- 

We can simplify this by adding the t times the pth column 
to this, yielding for the modified matrix 

M a ,p-1 = S a ,p-1 + t5 ap + 2iG a .2p-2- 



Continuing in this fashion gives 
p 

M ab ^Y. 6 *^ + 2lG «>»+P-^ 
3=b 

for b = 2 ... p and a = 1 . . . 2p. We already have 

Mai = iG a .2p- 

We can do the analogous manipulations for the columns 
c = p + 1 . . . 2p — 1. This yields 

c 

M ac = ^ S aj t c a — 2iG a , c -p+i- 
j-p+i 

and 

M a .2 P = -iG aA . 

We can now write M in terms of the Q s and R s defined 
in the last section. For example, 

fi(-l)[°/ 2 ]Qi_„ l<a<p 
a - 2v \(-l)^-P-^R a _ p p+l<a<2p 

where [x] is the greatest integer less than x. Putting this 
all together yields, for example, for p = 4 



M = 



— i?4 


2R 1 


2R 2 


-2R 3 


2iQ 1 





-2iQ 3 





-Rs 


1-2R 


2i?i 


2R 2 





2iQ 1 





-iQi 




t-2R_! 


l-2R 


2i?i 


-2iQ x 





2iQi 





Ri 


t 2 + 2R_ 2 


t-2R- X 


l-2i? 





-2tQi 





iQ 3 







-2iQ x 





1-2J2q 


t-2i?_i 


t 2 + 2i?_ 2 


Ri 





2iQi 





-2iQi 


2Ri 


l-2Ro 


t-2R_! 




-iQi 





2iQ 1 





2R 2 


2Ri 


l-2i? 


-R 3 





-2iQ 3 





2zQi 


-2R 3 


2R 2 


2R X 


— i?4 



where we have used the facts that Q_ s = Q s and Q s = for s even. This matrix looks much nicer if we reshuffle rows 
and columns. We permute the first column through the others so that it becomes the pth column, and permute the 
2pth column through so that it becomes the p+ 1th column. This does not change the determinant. We then relabel 
the indices p + 1 . . . 2p in reverse order (i.e. interchange 2p <-> p + 1, 2p — 1 <-» p + 2, etc.) We also multiply the pth 
and the p + 1th columns by 2, so that we need to divide the resulting determinant by 4. These manipulations put the 
matrix into the form 



(Bl) 
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where 1Z and Q are the p x p matrices defined by 

Kij = (-1)^-^2^.^+9(1- j)^- 1 (B2) 
Qij = t(-l)^^2Q p+1 ^ 3 (B3) 

where 0(x) — 1 for x > and for x < 0. For example, for p = 6 we have 



ft = 



/ 


2i?! 




2R 2 


-2R 3 


-2R 4 


2R 5 


2i? 6 \ 


1 


-2i? 




2i?i 


2i? 2 


-2R 3 


-2i? 4 


2R 5 


i- 


- 2i? ! 


1 


-2R 


2R 1 


2R 2 


-2R 3 


~2Ri 




+ 2i?_ 2 


t- 


- 2i? ! 


1-2R 


2R l 


2R 2 


-2R 3 




+ 2i?_ 3 


t 2 


+ 2i?_ 2 


t - 2i?_i 


\-2Rq 


2R 1 


2R 2 


V 4 


- 2i?_ 4 


t 3 


-2i?_ 3 


t 2 + 2R- 2 


t - 2i?_i 


l-2R 


2Ri J 



Q = 



f 2iQ 5 





-2iQ 3 





2iQ x 


o \ 





-2iQ 3 





2iQ 1 





-2iQ l 


-2zQ 3 





2iQ 1 





-2iQi 








2iQi 





-2iQi 





2iQ 3 


2iQi 





-2iQ 1 





2iQ 3 





V o 


-2iQ 1 





2iQ 3 





-2iQj 



The determinant of matrices of the form (jBlfl can be simplified^ 
det^ 2 Vdct^ +S S Udct^+ 2 



= det(ft+Q) det(ft-Q). 



In our case, the latter two determinants are the same, 
because 1Z is symmetric around the off-diagonal (the di- 
agonal from upper right to lower left), while Q is anti- 
symmetric around the off-diagonal. Thus 

detM = i (det(K+ Q)f . 



This yields the result ll'.'il) for the monomer- monomer cor- 
relator above. 
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